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Abstract. Recent data, e.g. from SDSS and 2dF, reveal a robust bi-modality in the 
distribution of galaxy properties, with a characteristic transition scale at stellar mass 
M, ^ 3 X 10^" (near L,), corresponding to virial velocity V ^ lOOkms"^. Smaller 
galaxies tend to be blue disks of young populations. They define a "fundamental line" 
of decreasing surface brightness, metallicity and velocity with decreasing M*, which 
extends to the smallest dwarf galaxies. Galaxies above the critical scale are dominated 
by red spheroids of old populations, with roughly constant high surface brightens and 
metallicity, and they tend to host AGNs. A minimum in the virial M/L is obtained at 
the same magic scale. This bi-modality can be the combined imprint of several different 
physical processes. On smaller scales, disks are built by cold flows, and supernova 
feedback is effective in regulating star formation. On larger scales, the infalling gas 
is heated by a vmal shock and star formation can be suppressed by AGN feedback. 
Another feedback mechanism - gas evaporation due to photo-ionization - may explain 
the existence of totally dark halos below V ~ 30kms~^. The standard cooling barriers 
are responsible for the loose upper and lower bounds for galaxies: IQ <V < 300 kms~^. 

1 Introduction 

Observations indicate four characteristic scales in the global properties of galax- 
ies, expressed either in terms of virial velocity V , or stellar mass Af,: 

• An upper limit at ^ 300kms~^, which is the lower bound for clusters. 

• A robust bi-modality scale at Af, ~ 3 x 10^" Af© corresponding ioV^ 
100 kms^^. This is comparable to L» of the Schechter luminosity function, 
marking the exponential upper bound for disk galaxies. Smaller galaxies tend 
to be blue, star-forming disks with the surface brightness and metallicity 
decreasing with decreasing A^f, (LSB). More massive galaxies are dominated 
by red, old-population spheroids with high surface brightness and metallicity 
insensitive to Af* (HSB) ^3]. AGNs typically reside in halos above the bi- 
modality scale 16 . The halo M/L shows a minimum at the same scale |35| . 

• Dark-dark halos (DDH) with no luminous components seem to dominate 
the halo distribution below V ^ 30kms~^. The relatively flat faint-end 
luminosity function and the steeper mass function predicted by the standard 
ylCDM cosmology cannot be reconciled with the virial theorem and the 
TuUy-Fisher relation obeyed by the luminous galaxies {^^. Also, most dwarf 
spheroidals in the Local Group lie in this range, while dwarf irregulars are 
typically V > SOkms-^ ^ [TS| El- 

• A lower bound at ~ lOkms"'^ for dwarf galaxies (see [H]). 
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We summarize here how theory predicts the following characteristic scales: 

• A loose upper bound for cooling (by bremsstrahlung) on a dynamical time 
at M ~ lO^^-i^M©, corresponding to 300 kms"^ 

• A virial shock-heating scale at M* ~ 3 x 10^^ Mq, or V ^ lOOkms"^. In less 
massive halos virial shocks cannot form and the gas flows cold into the disk 
vicinity (0. 

• A supernova-feedback scale a.t V ^ lOOkms"^, below which the potential 
wells are shallow enough for the energy fed to the ISM by supernova to 
significantly heat the halo gas, or even blow it out i^^- 

• A photo-ionization scale at V ^ 30kms~^. Once the gas is continuously 
ionized by the background UV flux, IGM gas cannot accrete into smaller 
halos and gas already in such halos evaporates via a steady thermal wind 

• A lower bound for cooling (by atomic hydrogen) at y lOkms ^. 

We briefly address the origin of these scales, and try to associate them with 
the observed scales. 

2 Shock-Heating Scale 

The standard picture of disk formation is that after the dark halo virializes, the 
infalling gas is shock heated near the virial radius. It then cools radiatively and, 
if the cooling time is shorter than the dynamical time, the gas accretes onto a 
central disk and forms stars with the associated feedback. The upper-limit halo 
mass for efficient cooling is Mcooi ^ 10^^"^'^ Mq [321 . The common wisdom 
used to be that this explains the upper bound for disk galaxies, but it has been 
realized that the predicted scale is somewhat higher than the observed scale. 

Fig. n shows the time evolution of the radii of Lagrangian gas shells in a 
spherical system consisting of gas (with primordial composition) and dark mat- 
ter, as simulated by Birnboim & Dekel |2j using a ID hydro code. Not shown are 
the dark-matter shells, which collapse, oscillate and tend to increase the gravi- 
tational attraction exerted on the gas shells. The left panel focuses on massive 
halos of ^ Indeed, a shock exists near the virial radius at all times; it 

gradually propagates outward, encompassing more mass in time. The hot post- 
shock gas is in a quasi-static equilibrium, pressure supported against gravity. 
The right panel focuses on halo masses smaller by an order of magnitude. An 
interesting new phenomenon is seen here, where a stable shock forms and prop- 
agates from the disk toward the virial radius only after a mass of ~ lO^^Af© has 
virialized. In less massive systems, the cooling rate is faster than the infall time, 
making the "post-shock" gas unstable against gravitational infall and thus not 
being able to support an extended shock. 

This behavior can be understood via the following instability analysis of 
the gas just behind the shock 5 . The standard equation of state is written as 
P = (7 — l)p e where p and e are the gas density and internal energy, and 7 = 5/3 
for a mono-atomic gas. In the no-cooling case, the adiabatic index is defined as 
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Fig. 1. Time evolution of the radii of Lagrangian gas shells (lines) in a spherical sim- 
ulation of a protogalaxy consisting of primordial gas and dark matter. Temperature is 
marked by color. A shock shows up as a sharp break in the flow lines and as an abrupt 
change in temperature. The lower discontinuity marks the "disk" radius, formed due 
to an artificial centrifugal force, (a) A massive system, where the virialized masses are 
in the range lO" - lO^M©. (a) A less massive system, 10^'' - 10^^ Af©. A virial shock 
exists only in systems more massive than > W^^Mq. 



7 = (91nP/91np)ad, and the system is known to be gravitationally stable once 
7 > 4/3. When there is radiative cooling at a rate q, we define a new quantity: 

d\nP\ _ pq _ tdyn 



ydlupj pe tcooi 

The second equality follows from energy conservation, e = —PV — q, plugged 
into the equation of state. Note that 7cff = 7 when g = 0, and that the difference 
between the two is nothing but the ratio of characteristic times idyn/^cooi, where 
tdyn = p/p — fs/u, with Ts the shock radius and u the infall velocity into it. 
Using the jump conditions across the shock, one can express 7cff in terms of the 
pre-shock gas quantities: 

(7 + 1)^ pr.A{T) ^_ iJi 27-1 , 
^''^"^ 6(7-1)2 Ki ' ' -fcBiVA(7 + l)2" ' ^ ' 

where p is the pre-shock density. 

In order to test for stability, we perturb a shell hy r ^ r + 5r and compute 
the sign of the restoring force: 5r / Sr. We assume a homologous infall 5r = u5t 
(supported by the finding in the simulations above), and balancing forces f = 
-p-^VP - GM/r^ = 0. This leads to 

5r/5r cx 7 - 4/3 - 27^-^^1(7 - 7eff) . (3) 

The system is stable when this expression is positive, namely for 



(4) 
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For 7 = 5/3, we get 7crit = 10/7 replacing the 7crit = 4/3 of the adiabatic 
case. The critical scale is then characterized by icooi/idyn = 21/5 at the disk 
vicinity, which translates to icooi/^dyn — 1 at the virial radius. Note that in our 
notation tdyn — R/V at the virial radius, while the dynamical time used in the 
order-of-magnitude arguments of [2Z| and is several times larger. 

When 7efi is computed using eq. Q just above the disk radius at every time 
of the spherical simulation, we find that the criterion of eq. Q indeed works very 
well. Before the shock forms, it is well below 7crit and gradually rising, reaching 
7crit almost exactly when the shock starts propagating outward. It then oscillates 
about 7crit with a decreasing amplitude, following the oscillations in the shock 
radius seen in Fig. Q 

When put in a cosmological setting, we find that the critical halo mass for 
virial-shock heating is 

M.hock ^ (1 - 6) X 1O"M0 . (5) 

The higher estimate is obtained for near solar metallicity. This scale is compa- 
rable to the observed bi-modality scale, and is quite insensitive to redshift. To 
see this, we use the spherical-collapse virial relations at z: M oc (1 -I- z)~^/'^V^ oc 
(1 -I- z)^R^, and the general trend in the relevant range A oc T^^/^, to obtain 
icooi/^dyn OC M independent of z. Since at 2 > 2 most halos are of M < Mshock, 
we conclude that in most forming disks the gas never heats to the virial temper- 
ature - it flows cold all the way to the disk vicinity. 

Although we are back to a condition involving ^cooi/^dyn, which seems to re- 
mind us of the original criterion by Rees & Ostriker, it is now based on analyzing 
the actual physical process that could be responsible for the heating to the virial 
temperature, namely the existance of a virial shock. Because the dynamical time 
involved in the current analysis is smaller than the one used in the original es- 
timate by a factor of a few, our derived upper bound for halo masses in which 
cold flows prevail is smaller by a similar factor. 

Preliminary results from cosmological hydro simulations indicate that the 
phenomenon is not restricted to spherical systems, and that our stability criterion 
may in fact be approximately valid in the general case where the inflow is along 
filaments. Fig.[51shows snapshots from Eulerian simulations by A. Kravtsov 19.^ 
|2()| . showing the gas temperature and flow fields in two protogalaxies: one of 
M ~ 3 X 10"Mq at z ~ 4, and the other of A/ ~ 2 x 1O^°M0 at z ~ 9. While 
the more massive galaxy shows a hot gas component at the virial temperature 
behind a virial shock, the smaller galaxy has mostly cold gas inside the virial 
radius. Also seen in the massive galaxy are cold streams penetrating the hot 
medium toward the central disk. Similar results have been obtained by Fradal et 
al. , who emphasis the feeding of galaxies by cold flows preferentially at early 
epochs, corresponding to less massive galaxies. We are in a process of testing 
our shock-stability criterion using cosmological simulations, and generalizing it 
to the realistic case of elongated streams. 

Beyond introducing a new scale in galaxy formation, the cold streaming to- 
ward the disk, as opposed to the commonly assumed shock heating and slow 
accretion, may have several interesting implications. For example: 
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Fig. 2. Snapshots from cosmological hydro simulations (courtesy of Audrey Kravtsov), 
showing the gas temperature and flow fields in two protogalaxies. (a) A halo of Af ~ 
3 X 10^^ Mq at z ~ 4, showing a virial shock and heated gas penetrated by cold streams, 
(b) A halo of M ~ 2 X lO^^M© at z ~ 9, showing cold infall with little shock heating 
inside the virial radius of ~ 7/i~^kpc. 



• It may modify the star formation rate and the associated supernova feedback. 
If the result of the cold stream hitting the disk is a burst of stars at high z, 
this may help the semi-analytic models (SAMs) explain better the presence of 
luminous red galaxies at z ~ 1. This process should be properly incorporated 
in the SAMs. 

• It may reduce the soft X-ray flux from halos. This is consistent with observed 
individual galaxies 6 , and may explain why the observed X-ray background 
is lower than the predictions based on shock heating in all halos |26j . 

• The cold gas may emit a significant fraction of its infall energy in Ly-a, 
either if the streams cool and merge into the rotating disk gradually without 
shocking iITj, or after shocking at the disk vicinity. In the latter, the X-rays 
produced are likely to be captured inside a dense Strongrem sphere of a few 
kpc and the energy may eventually be released by Ly-a. This may explain 
the observed Ly-a emitters (21] . 

3 Supernova Feedback Scale 

Supernova feedback, which is a very different physical process, gives rise to a 
characteristic scale in the same ballpark. The energy fed to the ISM by super- 
novae can be written as 

Es-N ^ l^eMJrad OC M,(tcoolAdyn) , (6) 

where v is the number of supernovae per unit mass of forming stars, e is the typ- 
ical supernova energy, and icooi is the time available for the supernova remnant 
to share its energy with the medium during its adiabatic phase, before a signifi- 
cant fraction of it is radiated away. The second proportionality follows from the 
crude assumption Af» ^ M^/t^yy^. Dekel & Silk |S] pointed out that for halos in 
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the relevant range, where the virial temperature is T lO^K, the cooHng rate 
behaves roughly like A oc T~^, and then tcooi/idyn — 0.01 for all galaxies. This 
means that, despite the significant radiative losses, the energy fed to the ISM 
is i?sN oc M,. When this energy input is compared to the energy required for 
significantly heating the gas or blowing it out, -^binding — Mg^.^V'^, one ends up 
with the critical scale for effective feedback, 

VsN =i lOOkms-^ (7) 

In shallower potential wells, the feedback from a burst of stars can significantly 
suppress further star formation, and produce low-surface-brightness (LSB) and 
dwarf galaxies. We note that Vsn corresponds to Af, ~ 3 x lO^^M©, comparable 
to the observed bi-modality scale. 

Tight correlations between the properties of galaxies below 3 x IO^^Mq define 
a "fundamental line" (FL) which stretches across five decades in luminosity 
down to the smallest dwarfs. The mean scaling relations, involving M^,, V, the 
metallicity Z and the surface brightness ^, are approximately jJS] [S] [12 |U 

V^M°-^, Z<^M°-\ (8) 

Note that the familiar TuUy-Fisher relation at the bright end, V oc M^-^ 0, 
seems to fiatten in the dwarf regime. Supernova feedback can explain the origin 
of the FL in simple terms [H]- The above energy criterion, E'en Af, oc Mg^^gV^, 
with the assumption that the relevant gas mass is a constant fraction of the halo 
mass, Mgas oc M, immediately implies the relation driving the FL: 

M,/M oc . (9) 

Assuming spherical collapse to virial equilibrium, the halo virial quantities are 
related via, M oc V'^ oc R^. In the instantaneous recycling approximation, we 
approximate Z oc M,/Mgas- For the stellar radius we adopt the standard assump- 
tion of angular-momentum conservation |24j . oc XR, with A a constant 
spin parameter, and use it in /.t oc M^j R^. Combining the relations from eq. ^ 
and on, we recover the scaling relations, eq. ©. This success of this simplest pos- 
sible toy model indicates that supernova feedback could indeed be the primary 
driver of the FL below the bi-modality scale, and it strengths the association of 
this process with the origin of the transition scale itself. 

4 The Origin of Bi-Modality 

In Fig.m we schematically put together the physical scales discussed so far in the 
classic density-temperature diagram. Galaxies are loosely confined to the region 
where tcooi < iff - The line M^^shock ~ 3 x lO-'^^M© distinguishes between massive 
galaxies which suffer virial shock heating and lower mass galaxies where the gas 
flows cold toward the disk. The line V ^ lOOkms"^ marks the transition from 
deep potential wells where HSBs form to shallow potential wells where supernova 
feedback is effective and leads to LSBs and dwarfs along the FL. 
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Fig. 3. the characteristic scales in the classic diagram of gas density versus virial 
velocity. The horizontal (gray) curve marks 2-a halos in /ICDM. The V-shaped (green) 
curve refers to tcooi = tdy-n\ galaxies can form only above this line, (a) The shock- 
heating scale, (b) The supernova feedback scale. The two scales are comparable, at 
V ~ 100 kms"^ or Af* ~ 3 x 10^° Mg. 



The two different physical processes happen to yield a similar characteristic 
scale, and they thus combine to produce the bi-modality. A possible scenario is 
as follows: 

• Halos below Mghock had cold infall to disks, possibly associated with bursts 
of star formation. Galaxies above this threshold are built from progenitors 
that were below the threshold at high z and thus formed early disks, which 
later merged to spheroids of red, old stars. Any halo gas in such massive 
halos consists of a hot medium and possibly cold streams. 

• In halos below Vsn, the SN feedback regulates the star- formation rate, which 
leads to blue, young stars and the FL, M^jM cx y^, with MjL rising toward 
smaller halos. The SN feedback may also suppress any AGN activity in these 
halos. 

• In more massive halos, the gas, which has been heated by virial shocks, may 
be prevented from cooling by AGN feedback, causing M/L to rise toward 
larger halos in proportion to the total-to-cold gas ratio. There seems to 
be enough energy in AGNs, though the actual feedback mechanism is not 
properly understood yet. 

• Protogalaxies below Afghock may be seen as Ly-a emitters rather then as 
X-ray sources, while the hot gas in more massive halos still produces X-rays. 

5 Photo-evaporation Scale 

Some of the halos below V ^ 30kms^^ must be totally dark, with no luminous 
component. Marginal evidence is starting to emerge from the statistics of image 
multiplicities in gravitational lenses (18) . and independently from the following 
comparison between theory and the observed frequency of dwarf galaxies. While 
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the faint-end galaxy luminosity function is roughly 4>{L) oc L^^ (or steeper), the 
halo mass function predicted in ylCDM is roughly ip{M) oc (or flatter). 

These cannot be reconciled with the virial relation M oc and the TuUy- 
Fisher relation in the dwarf regime, L (x (or flatter). If, however, only a 
fraction fhiV) of the small halos have a luminous component (with L/M oc 
as discussed above) , then the observations can be reconciled once /l oc ^ ^ . Thus, 
if the dark-dark halos appear at ~ 30kms~^, we expect ~ 90% of the halos 
to be totally dark at the faint end, V ~ lOkms"^. 

Supernova outflows alone are not likely to produce DDHs, both because they 
must leave a trace of stars behind, and because they tend not to be equally 
effective in all directions . Ram pressure due to outflows from nearby galaxies 
may remove gas in some cases but probably not enough for the purpose. 
On the other hand, the origin of DDHs could be explained by radiative feedback. 
The first stars and AGNs produce UV flux which ionizes most of the gas in the 
universe by Zion ~ 10, and keeps it at T ~ (1 — 2) x lO'^K until Zcnd ~ 2 or 
so Further infall of gas is suppressed in halos below the Jeans scale, which 
is crudely estimated to be about V ~ 30kms~^ but this does not explain 
the complete removal of the gas already residing in those halos. A significant 
fraction of T ~ lO'^K gas is expected to escape instantaneously from halos of a 
correspondingly shallower potential, V < lOkms"^ T, but not from halos of 
10<V < 30kms-^ 




Fig. 4. The evaporated mass fraction as a function of virial mass and velocity (insert). 
The curves on the left correspond to instantaneous expulsion, and the curves on the 
right refer to steady evaporation, where the ionization is assumed to be efficient from 
z = 8 till either 2 = 4 (solid) or 2: = 1 (dashed). The gas temperature is either lO^K 
(left) or 2 X lO^K (right). While half the gas is lost instantaneously from 1/ ~ 10 kms~^ 
halos, a steady wind can remove a similar fraction of the gas from V ~ 30kms~^ halos. 
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However, if a continuous input of energy keeps the gas ionized for many 
dynamical times, it can evaporate via a steady thermal wind even from deeper 
potential wells. Once the gas is lost from the shallow edge of the potential well, 
it is continuously replenished by interior gas, and so on. In Shaviv & Dekel [22! 
we analyzed this process following the classic solar- wind analysis ) and found 
that a significant fraction of the gas eventually evaporates from halos as large 
as y ~ 30kms~^. Fig. 0] summarizes the results of this analysis. 

Thus, halos in the range 10 < ^ < 30kms~i^ lose most of their gas after the 
cosmological reionization. They may end up as gas-poor dwarf spheroidals with 
an old stellar population ^ [231, or as DDHs (see also jSIIl)- 

Halos oi V < 10kms~"i are not expected to form stars at all because of the 
cooling barrier at T ~ lO'^K (Fig. O. Atomic line cooling is very efficient just 
above this temperature, but the only cooling agent below W^K is molecular hy- 
drogen, which cannot survive the UV background and is anyway very inefficient 



6 Conclusion 



The following table summarizes the characteristic scales originating from the 
physical processes discussed above and their possible association with the ob- 
served scales. The classic argument of cooling on a vaguely-defined dynamical 
time provides relatively loose upper and lower bounds for luminous galaxies. 
The combination of (a) the newly discovered scale distinguishing between virial 
shock heating and cold flows and (b) the familiar supernova feedback scale seems 
to produce the robust bi-modality imprinted on the observed galaxy properties, 
and in particular the characteristic upper limit for disk galaxies near i*. This 
is provided that AGN feedback, or another process such as thermal conductiv- 
ity, indeed prevents the shock-heated gas in massive halos from ever forming 
stars. The details of the proposed scenario are yet to be worked out (see also |3] 
^). While supernova feedback can be the primary process responsible for the 
"fundamental line" of LSB and dwarf galaxies, the radiative feedback due to cos- 
mological reionization may be responsible for totally evacuating small halos from 
gas and producing dwarf spheroidals and dark-dark halos. Once these physical 
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process are properly incorporated in semi-analytic models of galaxy formation, 
they may help solving the apparent conflicts between theory and observation. 

Acknowledgment: I thank my collaborators Yuval Birinboimi and Nir Shaviv, 
and those responsible for the 3D hydro simulations, R. Dave, N. Katz, D. Keres, 
A. Kravtsov, & D. Weinberg. This research has been supported by ISF 213/02 
and NASA ATP NAG5-8218. 



References 

1. Barkana R., Loeb A., 1999, ApJ, 523, 54 

2. Barkana R., Loeb A., 2001, PhR, 3 

3. Benson A. J., Bower R.G., Frenk C.S., Lacey C.G., Baugh CM., Cole S., ApJ, 599, 
38 

4. Binney J., 2004, MNRAS, |astro-ph/0308172| 

5. Birnboim, Y., Dekel, A., 2003, MNRAS, 345, 349 

6. Ciotti L., Pellegrini S., Renzini A., D'Ercole A., 1991 ApJ, 376, 380 

7. Courteau, S., MacArthur, L.A., Dekel, A. et al., 2004, astro-ph/0310440] 

8. Dekel, A., Silk, J., 1986, ApJ, 303, 39. 

9. Dekel, A., Woo, J., 2003, MNRAS, 344, 1131 

10. Fall M., Efstathiou G., 1980, MNRAS, xxx 

11. Fardal, M.A., Katz, N., Gardner, J. P., Hernquist, L., Weinberg, D.H., Dave, R., 
2001, ApJ, 562, 605 

12. Gnedin N.Y., 2000, ApJ, 542, 535 

13. Grebel E.K., Gallagher J.S., Harbeck D., 2003, AJ, 125, 1926 

14. Haiman Z., Rees M.J., Loeb A., 1996, ApJ, 467, 522 

15. Kauffmann et al. 2003, MNRAS, 341, 54 

16. Kauffmann et al. 2004, in preparation 

17. Keres D., Dave R., Katz, N., Weinberg, D.H. et al., 2004, in preparation 

18. Cohn J.D., Kochanek C.S., 2004, astro-ph/0306171 , 

19. Kravtsov A.V., 2003, ApJL, 590, LI ' 

20. Kravtsov A.V., Gnedin, O.Y., 2003, ApJ submitted fastro-ph/0305199|l 

21. Manning C., Stern D., Spinrad H., Bunker A.J., 2000, ApJ 537, 65 

22. Mateo M., 1998, ARA&A, 36, 435 

23. Mac low M.M., Ferrara A., 1999, ApJ, 513, 142 

24. Mo H.J., Mao S., White S.D.M., 1998, MNRAS, 295, 319 

25. Parker E.N. 1960, ApJ, 132, 821 

26. Pen U. L., 1999, ApJL 510, LI 

27. Rees M.J., Ostriker J.P, 1977, MNRAS, 179, 541 

28. Scannapieco E. Ferrara A., Broadhurst T., 2000, ApJL, 536, Lll 

29. Shaviv, N., Dekel, A., 2004, astro-ph/0305527 

30. Somerville, R.S. 2002, ApJL, 572, L23 

31. Tremonti et al. (SDSS), 2004, in preparation 

32. van den Bergh S., 2000, The Galaxies of the Local Group (Cambridge University 
Press) 

33. White, S.D.M., Rees, M.J., 1978, MNRAS, 183, 341. 

34. Woo, J., Dekel, A., Courteau, S., 2004, in preparation 

35. Yang X., Mo H.J., van den Bosch E.G., 2003, MNRAS, 339, 1057 



